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ABSTRACT 



The objective of this study is to analyze the behavior of 
a typical sh t p structural member exposed to fire. 



Included in this work are chapters pertaining to the 
modelling of the fire process and the structural member to be 
analyzed. After determining the expected temperature exposure 
of a shipboard structural member, a heat transfer analysis is 
carried out and the resulting load is subsequently used in the 
failure analysis. 



The major thrust of this study is the analysis of the 
failure mechanism of the collapse of an I beam. Two methods 
of analysis developed by Kecman , and Hayduk and Wierzbicki, 
are summarized and applied to the case of a fire load. The 
recuced moments due to fire exposure are determined and 
finally recommendations are offered for the use of this 
information in ship structural design. 
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Chapter 1 



iNTRnniir.TTnM 



Traditionally, Marine Engineers have handled the e-f-Fects 
o-f -fire on shipboard structures through the use o-f damage 
control procedures. Ideally, a structure should be designed 
to withstand anticipated -fire loads. The objective o-f this 
work is to apply the methods o-f post critical analysis o-f 
thin-walled structures developed by Mechanical Engineers to 
analyzing -fire loads on board ocean engineering structures, 
spec i -f i cal 1 y ship structures. 

A ship structure is a highly redundant structure. I -f the 
load bearing capacity o-f one member is diminished, then the 
surrounding structural members carry more o-f the «ippl i sd 
load. Shipboard -fires are char ac ter i st i cal 1 y highly localized 
phenomenon which may e-ffect only a small percentage o-f the 
entire structure. By examining the e-f-fect o-f -fire on a single 
structural member, one will gain insight into the behavior o-f 
the entire structure. Ultimately, ship structures may be 
designed to withstand the e-f-fects o-f typical shipboard -fires. 

It is the intent o-f this study to examine the e-f-fect o-f a 
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fire on a typical shipboard structure, namely the I beam. The 

major thrust will be the analysis of the failure mechanism of 

the beam. Two methods will be introduced which analyze the 
post yield behavior of an I beam. Figure 1.1 shows 
graphically how a fire load diminishes the mechanical 

properties of a beam. 

The material yields upon reaching the maximum moment in 
Figure 1.1. Consider the following relation: 

Stiffness = Moment/Curvature. 

Beyond the yield point, stiffness is less than or equal to 

zero. Normally a structural member remains in the elastic 

range, and the moment remains less than M , the fully plastic 

o 

yield moment. However, through exposure to fire, the 

stiffness will decrease so rapidly that the moment will fall 
below a safe range. 

The ultimate objective of this study is to show how much 
the strength of an I beam is diminished when exposed to fire. 
To this end, the structure and the fire process will be 
modelled; and the effect of temperature on the mechanical 
properties of a selected material will be given. After these 
models have been developed, a heat transfer analysis will be 
performed. As a result of this analysis, it will be shown 
that bending collapse is the primary mechanism of structural 
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•failure. Next, the failure mechanism analysis methods will be 
fully developed. Finally, recommendations will be offered for 
applying the results of this work to the design of an entire 
ship structure. 



Chapter 2 



MODEL nF...SIRIICTIIRF 



This chapter will develop the motivation towards choosinq 
an I beam as a typical shipboard structural element. 
Additionally, a material will be selected; and the temperature 
dependent properties o-f that material will be given. 



2^1-Bjir.k.gc.nund 

The principal aspects o-f ship structural design most 
commonly employed by Marine Engineers involves the application 
o-f rules establ ished by Societies such as the American Bureau 
o-f Shipping. These rules are based upon the primary loads 
experienced by a ship during its lifetime. These loads 
include: still water bending moments and shearing forces; wave 
induced moments and shearing forces; impulse loads; thermal 
loads; and dynamic loads. 

These loads which are unique to ocean going vessels are 
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There-f ore , 



o-ften di-f-ficult to accurately predict. There-fore, ship 
structural design is largely a process involving guesses 
derived -from intuition or empirical data. The primary stress, 
strain and displacement behavior o-f the ship structure are due 
to displacements limited to those parallel to the unloaded 
midsur-face o-f the plate. Thus, tradi t i onal 1 y hull strength is 
determined through an analysis of the midships section. 
Figure 2.1 illustrates a typical midships section. 



Today, Ship Structural Engineers have adopted a model 
wh i ch is motivated by considering the loads transmitted to the 
structure. Loads transmitted to a stiffened plate < by 
slamming, etc.) are in turn transmitted to the stiffeners and 
from them to the point supports, and so on. Eventually a load 
transmitting grillage is derived. Next, boundary conditions 
will be used to define if the spans are treated as a series of 
isolated points or continuous beams. In most cases the 
boundary points will be those which can be described as simple 
or clamped supports. 

Figure 2.2 shows examples of typical members of a ship''s 
structure. For the purposes of this study, an I beam will be 
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analyzed as its behavior is considered represen tat i ve o-f the 
open beams in Figure 2.2. Figure 2.3 gives the dimensions o-f 
the beam which herea-fter will be used as the model to be 
anal yzed . 



The -first moment o-f area o-f the I beam will be required 
in the anal ys i s . 



zdA = t^(b - t^) (a 



t ) 
w 




< 2.1 ) 



A = Sec t i on Area 



t = Thickness of the flanges 
f 

t = Thickness of the web 



w 

a = Width of flanges 



b = Height of web 



2 

Assume t is negl igible. Hence, the first moment 
be approximated by: 



of area may 



zdA = abt^ 




( 2 . 2 ) 



2— 3_Sfi.lfi.c±lon_o£_lda±e.iii.al 

Several materials have been used throughout shipbuilding 
history starting with wood, an d p r ogr e ss ing to iron and 
steel . Today several grades of steel and even titanium are 
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us&d in the shipbuildinQ industry. One 0 "f the most c ommon 
grades of steels used is ABS Grade A. Material Scientists, 
however , have not tested ABS Grade A in order to determine the 
temperature dependent material properties. Ordinary 
structural steels such as ASTM-A36, have been widely used as a 
test specimen. Table 2.1 compares the properties of the two 
steel s . 



Carbon , max Y. 
Phosphorous, max X 
Su 1 pher , max V. 
Yield Stress, ksi 



Table 2.1 
Steel Proper ties 



ASTM-A36 

0.23 
0.04 
0.05 
3S .00 



ABS Grade-A 

0.26 

0.04 

0.04 

34.00 (min) 



assumption that an 
without sacrificing 



It is considered a conservative 
analysis for ASTM-A36 can be used 
accur acy , < r ef . 1 ) . 



One of the major contributions to the loss of strength is 
the temperature dependence of the mechanical properties of the 
materials under consideration. Several researchers have 
predicted the de ter i orat i on of steel structural elements at 
elevated temperatures. A search of the literature reveals 
good agreement among those whom have studied the temperature 
dependence of material properties. 
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Harmathy and Stanzak (.rsi, 2) have carried out several 



experiments in order to provide to Design Engineers, primarily 
to those concerned with -fire endurance of structural elements, 
the elevated temperature mechanical properties of common 
commercial steel products. The temperature dependent 
expression for Young's Modulus of Elasticity, E, <psi) for 
ASTM-A36 steel is given here: 

E = 30*10^ - 9.3T^ <2. 3) 

O 

T = Temperature of heat exposure, Fahrenheit 

Steel loses strength as it is heated. Figure 2.4 illustrates 
the temperature dependent tensile strength of ASTM-A36 steel. 

o o 

Notice that between 1000 and 1200 F the tensile strength 
falls below the usual design stress. 

For the purposes of this study, 1 inear dependence of 
temperature i s assumed . The foil ow ing equation will be used 
in this anal.ys i s . 

0^ = -20. 5T + 40,000.0 <2. 4) 
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Chapter 3 



tlDDEJ Qg_FTRF PRnrF.q.q 



A survey o-f the literature reveals several methods of 
char ac ter i z i ng a fire. Harmathy offers a numerical technique 
which is summarized in Section 3.1. Most Fire Engineers use a 
time tefTiperature curve to define the character i st i cs of a 
fire. The development of an acceptable time temperature curve 
is given in Section 3.2. 



3_..l_Cha£;az;±£i:l^U.jc^_ci^_a_Shl^ibj3aj::d_£lx;je 

The heat exposure on a shipboard structure resulting from 
a fire is a function of several identifiable and incidental 
variables, <ref. 3). The three primary factors contributing to 
the severity of fully developed fires is the flow rate of air 
into a burning compartment, (compartment ventilation); the 
total fire load; and the thermal charac t er i st i cs of the 
compartment boundaries. 
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Numerous 



theore t i cal 



studies 



have 



de term i ned 



time-temperature curves oi typical -fires based on these 
factors, <re-f 3). Even with simpli-fying assumptions, these 
■factors are only calculable with a computer. Harmathy 
introduces three easily calculable -fire severity parameters 
which can be assumed to be accurate predictions o-f the 
contributing -factors to a -fully developed fire. These 
parameters are defined as follows: 

1. The duration of the fully developed fire, t; 

2. The overall penetration flux, q, that is the heat flux 

absorbed by the compartment boundaries, averaged 

spatially over the boundary surfaces, and} 

3. The average temperature of the compartment gases, 

<average "fire temperature"), T , averaged spatially over 

g 

the period of full development. 

It is useful at this time to introduce the three periods 
of freely burning compartment fires. See Figure 3.1. The 
focus of this study is on the period of full development, 
since 80-90 percent of the fuel energy is released during this 
period. 



A compartment fire is a fire which is enclosed in a shell 
of ar b i t r ar y shape. Typically on board a ship, ventilation is 
secured to the compartment with the fire, and free 



22 



1 









OPb 



0 '= 






hGUVcS: "\ 



f 

! 

i 



1 

I 



I 

I 



II 



communication to the surrounding spaces is not allowed. Fire 
load is a numerical expression which is a function of the 
available fuel in the compartment. The units of specific fire 
load is pounds of wood per square foot, <ref. 1). 

To compute the fire load of a particular compartment, the 
products of the weights in pounds of the material, and the 
heat capac i ty are summed for all the combustible materials in 
the compartment. The resulting expression, which has units of 
BTU's must be divided by 8000 BTU/lb. This result will provide 
the same available energy as the equivalent weight of wood. 

Equation <3.1) gives the equivalent weight of wood in 
pounds divided by the deck area of the compartment. 

f = Ew.h./ha^ <3.1) 

11 w T 

1 

2 

F = Specific fire load, Ib/ft 

U1 = Weight of individual materials, lb 
i 

H = Heat capacity of individual materials, BTU/lb 
i 

H = Heat capacity of wood, BTU/lb 
w 

2 

A = Area of deck, ft 
F 

Table 3.1 gives typical fire load values for various 
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shipboard compartments, <re-f. 4). 



Tabl e 3.1 

Speci-fic Fire Loading o-f Shipboard Spaces 

lbs. /square -feet 

Control Spaces 

Ulhee 1 house/Char troom 1.5 

Fire Control Stations 1.5 

Escape Routes 

Corr i dors 1 . 5 

Stairway Enclosures 1.0 

Accommodation Spaces 
Staterooms 

Fire Resistant Furnishings 3.0 

Combustible Furnishings 5.0 

Publ i c Spaces 

Fire Resistant Furnishings: 

Lounges, Restaurants, etc. 3.0 

Ferry L'essels 1.5 

Combustible Furnishings 5.0 

Sanitary Spaces 0.4 

Ser V i ce Spaces 

Galleys 10.0 

Pantries <no heating appliances) 4.0 

Food Concession (-ferry vessels- 

no combustible stowage) 1.5 

Workshops 10.0 

Storerooms - Combustibles 10.0 

Storerooms - Cleaning gear only 3.0 

Ship''s Laundry 10.0 

Laundry - Priuate Use 1.5 

Main Machinery and Cargo Spaces 10.0 

Auxiliary Machinery Spaces 5.0 

Tanks and Coids 0-0 
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It has been determined by 



several researchers that the 



•flow rate o-f air into a burning compartment, (ventilation) is 
just as critical to the severity of the fire as the fire 
load. The determination of compartment ventilation depends on 
how the rate of entry of air compares with a critical value. 
The rate of entry of air into a compartment can be expressed 
as ; 

u = c^A h (3.2) 
a 1 w w 



U = Flow rate of air, 1 b/hr 
a 

2 

A = Area of window, ft 
w 

h = He i gh t of window, ft 
w 

5/2 

C = 240 Ib/ft -hr 
1 

The critical value is defined as; 

(3.3) 



C =28.1 Ib/ft -hr 
2 

v= Material constant which depends on the furnishings of 
the compartment. For typical living spaces, 
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0.55 < < 0.90. 



The air -flow -factor is expressed as -follows: 



C = 



(U ) 
a cr 



(3.4) 



Thus, a -fire is ventilation controlled i -f I < 1.0; and a 
-fire is -fuel sur-face controlled i -f I >= 1.0. 



The average rate o-f burning o-f combustibles, R, has been 
determined during the period o-f full development based upon 
empirical data. For ^ >= 1.0 

R = (3.5) 



2 

C = 4.57 Ib/ft -hr 
3 



For £<1.0 



R = 0.163U 

3 . 



(3.6) 



The mass flow rate of gaseous products out of the 
compar tmen t is: 



U = R -t- u (3.7) 

g a 
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It -foil cws that ; 



Q = R(0.9326h, + 0.0686H.) 

12 



<3. 8) 



where ; 



Q = average rate of heat evolution from combustible 
mater i al s 

H 5760 BTU/lb, heat of combustion of char 

1 

H 11480 BTU/lb, heat of combustion of volatile 

2 

decomposition products 

S = the fraction of heat released by the combustion of 
the volatile decomposition products inside the 
compar tmen t 

S is calculated based upon an estimation of the average 
f 1 ame length, 1 ; 

Z = <3. 9) 

1/3 

C = 1 . 1 1 f t 
4 

Z = C^U (3.10) 

5 a 

1/3 1/3 

C = 0.366 ft(h) /<lb) 

5 
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Then , 



if il = H 6 = 1 

c 

if £ > H 6 = (H 

c c 

where H is the height oi the compartment, 
c 

The duration of a fully developed fire, x, 
as foil ows ; 



for C = 1 



T = C^/4) 



(3.11) 

(3.12) 



calculated 



(3.13) 



2 

C = 0.205 ft (hr)/lb 
6 

for C < 1 T = c /({)^ (3.14) 

6 



The other two parameters may be obtained numerically -from 
the -following equation: 



<3 = J (P - 

t 



U c (CT - 

g g g 







4 

T ^)] 
o 



(3.15) 



where ; 

A = Total area for the compartment 
t 

Q = Average rate of heat evolution (BTU/hr) 

o 

C = Specific heat of compartment gases <BTU/lb- R) 

9 
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o 

T = Temperature o-f combustion gases, R 
9 

o 

T = Temperature at time = 0, R 
o 

“8 2 o 4 

<^ = Ste-fan-Bol tzman constant, 0.171*10 BTLI/-ft -hr-[ R] 



^ = Empirical -factor, 1.05, dimensionless 



T 

g 




(T 

o 



+ ^(2TIC/7T) 



l/2j 4^1/4 



<3. 1 6> 



where : 

^ = Empirical -factor, 0.9, dimensionless 

2 

« = Average thermal d i -f -f us i v i ty , ft /hr 

o 

k = Average thermal conductivity, BTU/ft-hr- R 

Figures 3.2 and 3.3 show the dependence of the overall 
penetration flux, and average temperature of gases, 

respectively, on the air flow factor. These figures were 

2 

determined empirically for a specific surface of 0.63 ft /lb. 
Curves are given for three values of the specific fire load; 

2 

3.1 , 6.1 , and 12.3 Ib/f t . 

Figure 3.2 shows clearly that the maximum penetration 
flux, and fully developed fire duration occurs at ^ = 1. The 
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speci-fic -fire loading appears to have little effect on the 
maximum flux, and the duration is usually quite short. It may 
also be concluded from Figure 3.3 that the hottest fires occur 
in ventilation controlled compartments. It is possible at 
this time to draw some conclusions about shipboard fires. 
Below decks fires which do not communicate freely with the 
atmosphere are ventilation controlled. These fires are more 
destructive as the duration is greater, the average 
temperature is greater, and the effective heat flux is greater 
than for an above decks fire. 

McCarthy demonstrated the usefulness of Harmathy's 
technique for predicting fire severity of a typical berthing 
compartment. Figures 3.2 and 3.3 can be used for a typical 
berthing compartment, and in accordance with Table 3.1, a 

2 

specific fire load of 3.0 Ib/ft was used. He assumed that 
the ventilation is secured, and that one water tight door is 
open for access to firefighters. 

2 

A = 15.0 ft 
w 

h = 10 ft 
w 

U = 11384.2 Ib/hr 
a 

2 

A = 150 ft 
F 

9 =0.75 
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9483.75 



[U 3 
a cr 

i = 1.20 

Hence, the fire is fuel surface controlled. Notice that 

as A increase the fire approaches ventilation controlled. 
F 

From Figures 3.2, and 3.3; 

3 

T = 1.15*10 sec = 19.2 min 

2 

q = 6344 BTU/hr-ft 
o 

T = 1256 F 
9 

It should be appreciated from this discussion that fire 
is indeed a complex process with several unquan t i f i ab 1 e 
variables. The term destructive potential was introduced in 
reference 1. The destructive potential of a fire may be 
considered a function of the maximum temperature obtained, the 
duration of the fire, the duration of time at maximum 
temperature, and the quantity of energy transmitted from the 
fire. Most Fire Engineers use a time temperature curve as a 
definition of heat exposure. For the purpose of this study, 
the use of an accepted time-temperature curve will be 
considered an adequate expression of the fire process. 
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Traditionally, -fire tests have been per-formed based upon 
prior established standard time temperature curves. The most 
widely used -fire test in this country is the ASTM-E119 which 
was developed by the American National Standards Institute, 
<re-f. 5). This method vjas established as a result oi the 
necessity to standardize the structural performance of 
building construe t i on . 

The basic procedure of any fire test is to place the 
element to be tested into a gas furnace in which the standard 
time temperature curve is duplicated. The gas temperature is 
closely monitored through the use of thermocoup 1 es . 
Additionally, deflection, axial thrust, and collapse behavior 
are measured. The Ell? time temperature curve initially rises 

o 

rapidly reaching 1000 F at 5 minutes into the test. Beyond 
this point, the temperature is increased more slowly so that 

o 

at the end of one hour the furnace temperature is 1700 F. The 

o 

maximum test temperature, 2300 F is reached at 5 hours and 
held constant thereafter. 

Another test method commonly used by Ship Structural 
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EnQineers is the International Convention -for the Safety of 
Life at Sea <SQLAS) A60 test. This test is used for testing 
non-s true tur a 1 bulkheads from berthing and office 
compartments. The specimen observed is not loaded but 
strained at the edges. Other non-struc tural bulkhead tests 
are SOLAS A30 , A15, AO, B15, BO, but the requirements for 
these tests are not as severe as for the ASO test, <ref. 1). 

For the purposes of this study, reference 6 supports the 

o 

use of a step change to 1700 F vice the gradually increasing 
temperatures of ASTM-Ell? or the SOLAS ASO curves. Figure 3.4 
gives the ASTM-Ell 9 curve, the SOLAS ASO curve, and the 
assumed time temperature relation to be used in this 
analysis. 
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Chapter 4 



OEAI- IgAfciSEES-AMaLYS I S 



The intention o-f this analysis is to determine the 

relation between temperature and the distance along the length 

o-f the beam. A beam exposed to a -fire load will be modelled 

as a -fin with a temperature gradient. The temperature at the 

exposed face is T ; the temperature at the end of the fin is 

b 

T , the ambient temperature. The geometry of the fin is 
a 

defined in Figure 4.1 



4—J CDaduc±lDa_Jdaa±_Inan.s£axi 



The following assumptions are used: 

- T = Temperature at x = 0 

b 

- T = Temperature at x = L 

a 

- The heat transfer is steady state. 
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- k = Uniform conductivity of the material, 

o 

BTU/Cin-sec- F3 

- h = Uniform surface coefficient of heat transfer, 

2 o 

BTU/Chr-f t - F3 . 

Additionally, this study will use the assumption that the 
rod is thin, which implies that the temperature in the thin 
direction, perpendicular to x, is uniform. Therefore, the 
surface temperature equals the centerline temperature, (ref. 
7). Applying an energy balance on an elemental section *x , 
yields a governing differential equation. In steady state, 
with no internal heat generation; 

Expanding the first term of the right hand side gives the 
following equation: 

9q^ 

% ~ %,+Ax " (4.2) 

~ ^x+Ax (4.3) 

Sc ■ *^x+Ax " "|^(-kAyAz-^)Ax (4.4) 

Sc - Sc+Ax =|-(k&g)Ax <4.5> 

Where S = Surface area = AyAz. Conduction heat transfer across 
A A is given in the following 



40 



equat i on : 



II 

•?i>^ 


<4. 6) 


a 




to ■ - ’'a’9 


<4. 7) 


= hPAx(T - T ) 


<4. 8) 



Where 



P = Perimeter = 2<-iy + ziz) 

T = Temperature o-f -fin at x 

Substituting <q.5) and <4. 8) into <4.1) yields: 



^(ks|^)Ax - hPAx(T - T ) = 0 
dx dx a 



<4. 9) 



Finally, the governing differential equation for 6 = CT - 



T ] is; 
a 



— -^ 9.0 



<4.10) 



Where S, k , and T are constant. 

a 



The differential equation for one dimensional steady 
state heat flow with one internal heat generation source is; 



d^e . - «i 



(4.11) 



dx k 

Note that the fin problem is similar to the single internal 
heat 
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generation source problem with a negative heat source; 



^^2-Solu±xon 

The method o-f solution will follow a method provided in 
reference 8. Assume the solution is of the form given in 
equation (4.13): 



The boundary conditions for a fin of finite length, L, as 
defined by Figure 4.1 are; 




-Bx 



(4.13) 



Where 



B = [hP/kS]^^^ 



(4.14) 



X 



0 ; 



e _ T _ T B 9 
b a b 



(4.15) 



@x = L; 6 = T - T = 0 



(4.16) 



a a 



Substituting (4.15) and (4.16) into (4.13) yields 




(4.17) 



BL 



-BL 



(4.19) 
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Sol V i ng -for C : 
1 



r - 

Cl - b 



-BL 



BL -BL 
e -e 



C -follows from (4. 17), Finally the solution is: 
2 



<4. r?> 



e gB(L-x) ^-B(L-x) 

~ iZ 

b e - e 



<4.20) 



At this point, 
beam in Figure 2.3, 



it is necessary to solve -for B. For the I 
the perimeter to surface area ratio is; 



P 2(b + 2a - t ) 
i = ^ 

bt - 2t t + 2t.a 
w f w f 

Figure 4.2 gives the relation between thermal 
temperature for plain carbon steel, (ref. 9). 

o o 



<4.21 ) 

conduc t i V i ty and 
Within the range 



of interest, 70 F <= T <= 2000 F, the curve is considered 

approx i mate 1 y constant. For the purposes of this analysis, 

o 

the thermal conductivity at 1000 F wi 1 1 be assumed. 



k = 0.0005 BTU/in sec °F 



<4.22) 



The condition of natural convection in a gas is assumed 
for this analysis, <ref. 7): 

h = 2.0 BTU/hr ft^ °F <4.23) 

Theref ore , 

B = (P/10. 83)^'^^ <4.24) 

where, P/S is given in equation <4.21) 
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The validity ot' thin -fin theory must now be checked. 



According to reference 7, the thin fin theory is valid if: 



k > 
hb 



3 



(4.25) 



Substituting equations (4.22) and (4.23) into (4.25), the 
maximum value for b is; 



b = 3.6 ft 



(4.26) 



This is a reasonable criteria for this study. Therefore thin 
fin theory is valid. 



4^3_fiasu±±s 

It is now possible to draw some conclusions about the 
heat transfer through the beam. A Microsoft BASIC computer 
program was written in order to examine several geometric 
configurations of the beam. The source code for the program 
is in Appendix A. The results are given graphically in Figure 
4.3. The dimensions of the beams examined are in Figure 4.4. 

.The following conclusion may be drawn from Figure 4.3. As 
P/S, and L increase, the effect of the fire becomes more 
localized. With a typical shipboard I beam, the ratio of 
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perimeter to surface area is such that the fire load may be 
treated as a point load. That is to say that the heat 
transfer through an I beam remains highly localized in the 
vicinity of the fire. All effects of the fire virtually 
disappear at a distance which is a quarter of the length of 
the beam away from the fire. 
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Chapter 5 



AM^^l Y.qig np THF FAT I IIRF MFPHANTgM 



As was stated in the preceding chapter, an I beam exposed 
to a typical shipboard -fire will exper i ence the same -failure 
mechanisms as a beam under a point load. Recently there has 
been increasing interest in this type o-f -failure due to 
mechanical loads. It is the intention o-f this study to 
approximate the shipboard -fire load as an additional 
mechanical load, and to use two o-f the recently developed 
techniques -for analyzing the failure mechanism. Additionally, 
this chapter will include summaries of literature which 
focuses upon bending collapse of thin walled structures. 

Kecman <ref. 10) studied the bending collapse behavior of 
rectangular and square section tubes both theore t i cal 1 y and 
exper i men tal 1 y . He used a limit analysis technique and a set 
of formulae to derive the relation between the hinge moment 
and the associated angle of rotation. 'v'ery good agreement was 
found between his theoretical predictions and the experimental 
results. His theoretical predictions were in the form of a 
numer i cal solution. However, Hayduk and Wierzbicki (ref. 11) j 
gave a closed form analytical solution. In reference 11, 
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lower end upper bound solutions "for the mean crushing strength 
o-f cruc i -forms were obtained by considering bending and 
extensional modes o-f de-formation. 



Kecman developed a numerical scheme to analyze the 
bending collapse o-f rectangular and square section tubes, 
<re-f. 10). The primary failure mode is local buckling and 
hinge development. He considered only uniaxial bending 
col lapse. 



5.1.1 Maximum Bending Strength 

Kecman referred to a study which examined the bending 
strength of a box shaped section with the webs in tension and 
the top flange in compression. The effect of the webs was 
found to be negligible. Therefore, the critical stress of an 
I beam may be assumed to be the same as the critical stress of 
a box shape beam. 

The critical stress, a , of the section was derived in a 

cr 
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form : 



a 

cr 



2 

7T E 



12(1-V ) 



■2-(5.23 + 0.16^) (-) 

D 3. 



<5. 1 ) 



a = -flange width 



b = web height 



t = web, and flange thickness 



If the critical stress is less than the yield stress, a of 

o 



the material , then the flange buckles. The effective width 
becomes : 



a 

a = a(0.7-^ + 0.3) 
e a 

e 



<5. 2) 



cr = stress at the edge of the flange 
e 

If the critical stress is greater than or equal to the yield 

stress, then a = a. For the purposes of this study, assume 

e 

that collapsing occurs when a = a . It follows that the 

e o 

maximum moment, M , may be calculated using (5.1) and (5.2). 

m 



I f CT <= a : 
cr o 



M = a tb *- 
m o 



2a + h + a (3^+2) 
e b 



<5. 3a) 



3(a + b) 



If <r >= 3 o- ; 

cr o 



Mjj, = Mq " - t) + 0.5(b - 2t)^] 



<5. 3b) 
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a < o 



< 3 (7 , a mofTie n t , M 
o 



/ 



is defined; 



o 



cr 



o 



M ' = a tb(a + ^) 
o o 3 



<5. 3c) 



With the resulting maximum moment 



= M ’ + (M 
m o ' ( 




(5.3d) 



o 



5.1.2 Actual Bending Collapse Mechanism 

A detailed analysis of the hinge deformation pattern 
offers insight into the failure mechanism. The deformed 
pattern is given in Figure 5.1. 



mechanism were given in reference 10; 

(a) Most of the plastic deformation is concentrated along 
the stationary yield lines; JF , HI, EB, GB, FC , HC , CR, 
BA, GK, and EL. 

<b) Travelling yield lines of rolling deformation were 
also observed along GA , AE, KA, and LA while the distance 
BA changes during hinge rotation. 

<c) The corner angle along KL remains almost constant. 



The following characteristics of the deformation 
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<d) The de-formation o-f the tension -flange includes the 
travelling yield lines NP , and MQ, approx i mate 1 y 
corr espond i ng to the intersection o-f the unde-formed 
■flange planes and the current cylindrical surface TU. See 
Figures 5.1, and 5.2. 

<e> In-plane deformation of the Malls was also observed, 
particularly along the rolling lines GA, AE, BA, HC, CF, 
IR, and JR. This deformation is essential to maintain 
hinge k i nemat i cs . 

<f) The rolling radius, r, varies slightly along GA , 
measured by cutting at right angles to GA and its 
estimated value was approx imate 1 y proportional to the 
h i nge 1 ength . 



r =~ <0.03-0.05)*KL 

The radius reduces during collapse, so that rolling 
becomes more and more difficult, until the mechanism jams 
and starts a new phase of collapse. 

<g) The rolling radius along KA, and LA varies from 
infinity at K and L to a value close to r in (f) above. 



Four distinct phases in the hinge development were 
i den t i f i ed ; 

<1) Initial phase with a bulge in the web, A , in Figure 

1 
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5.3. 



<2) Second phase with rolling as described by the 

cross-section shape A in Figure 5.3. 

2 

(3) Third phase indicated by a jamming of the rolling and 

described by the cross section B A D in Figure5.3. 

3 3 

(4) Final phase initiated by the contact of the two 
buckled halves of the compression flange. There is no 
clear transition from one phase to another. 



5.1.3 Theoretical Bending Collapse Mechanism 

Since structures most commonly experience hinge collapse 
described by the second phase which includes rotation angles 
in the range, 5-10 <= <= 25-35. Kecman based his theoretical 

analysis on the following assumptions: 

- The web and the flanges of the section deform along the 
concentrated yield lines only; 

- The web and flanges are incompressible and 
i nex tensi bl e ; 

- The structural continuity is maintained in two 
charac ter i st i c sections; and 
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- Some in plane de-formation is essential -for hinge 
k i nemat i cs. 

The model is given in Figure 5.3. 



The rolling radius is negligible in comparison with the 
other hinge dimensions. Points A, B, C, and R are in a 
straight line. The coordinates o-f point B are: 



B(x,y,z) = (h, bcosp - [bsinp(2h- bsinp)]^'^^, 0) <5. 4) 

From continuity of the middle section, the web height, b, is 
determined: 



/ 2 2 , 1/2 



<5. 5) 



With y = y 5 and equations <5. 4), and <5. 5): 
A B 



z = bsin^p - hsinp + [bsinp(2h - bsinp) ] *cosp 
A 



An analysis of the rotation of the fiber in the top 

flange will lead to the equation for the hinge length, 2h . For 

a fiber which is initially parallel to the beam and passes 

through point A, longitudinal continuity is implied. The 

distance from this fiber and the edge KL is AD = b - z . 

A 

Initially, the coordinates of the fiber and the line GK are 

CO, b-z , 03. After rotation, the coordinates are [<b 
A 

z )sin/7, <b - z )cosp, 03. Point A has the coordinates Ch, y , 
A A A 
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z ]. The distance from point A and the rotated point is h if 
A 

the web and flanges are i nex t ens i b 1 e . The equation for h is 

ther ef ore ; 

= [h-(b-z^)sinp]^ + [y^- (b-z^) cosp] ^ ^ A (5.7) 

Substitution of <5. 4), and (5.6) into <5. 7); and meeting the 
condition required by <5. 5), it may be said for any ratio a/b, 
and any rotation angle: 

2h = a or 2h = b <5.S) 



The value which yields the lowest energy absorption will be 
used. Thus, 

< 

2h = a if a = b 

<5. 9) 

2h = b if b = a 

The second phase of deformation is complete when the 

rotation angle equals the Jamming angle, a . Jamming starts 

J 

when the two buckled halves of the compression flange contact 
each other. From Figure 5.3 the jamming angle is; 



^ • /h- 0.5t^ 

Ct = 2arcsinp( — r ) 

j b 



<5.10) 



5.1.4 Nominal Energy Absorbed During Hinge Rotation 



The assumptions listed in the previous subsection imply 
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that all the energy is plastic and is absorbed along the yield 
lines o-f the theoretical collapse mechanisms. Relative 
rotation along the -flange edges JF, and IH results -from hinge 
rotation, see Figure 5.1. From geometry, o-f Figure 5.3: 

7T b 

Y = -j - P - arcsind - |^inp) <5.11) 



The energy absorbed along these edges becomes: 



W = W + W = 2ra a[-^ - p - arcsind - ^inp) ] <5.12) 

X ur XH p z n 

m = Average -fully plastic moment o-f a unit width o-f 
P 



the -f 1 ange 



Similarly, it may be said: 

W_ = VJ = m a[TT - 2arcsin(l - ^inp) ] <5.13) 

^ wrv p il 

The energy absorbed by the yield line AB can be -found 

-from <5.11), where the length o-f AB = z is known -from <5. 6), 

A 

and the relative rotation is <. ^ - 2^). Hence: 

2 1/2 b 

= W = m [bsin p-hsinp+ (bsinO (2h-bsinp) ) *cosp ] [TT-2arcsin (l-r^inp) ] 
j AB p h 

Similarly, along BG, the length is h, and the relative 
rotation is ^/2i 



w 



4 



W. 



BG 



+ W. 



BE 



m hTT 

P 



(5.15) 



Figure 5.4 shows the relative rotation, ^/ , along yield 
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line GK . The construction o-f this -figure is detailed as 



follows: the plane co is normal to GK and contains point A. 

xy 

Point A'' is the projection of A onto the xy plane. A is 

o 

obtained by rotating abou t its intersection w with the xy 

xy 



plane. Coordinates of A" are: 



= 



htanp + bcosp - [bsinp (2h-bsino) ] 

, 2 
1 + tan p 



1/2 



<5.16) 



^A" " 



<5.17) 



The energy absorbed along both the GK, and EL 1 ines are then: 



W = + W„ = 2m barctan [• 

b GK EL p 






<5.18) 



The rolling radius along the yield lines GA, and AE 
remains relatively constant, so: 



W — W + w = 2m ( — )z 
6 GA AE p r A 



<5.19) 



During hinge rotation, the rolling radius varies as 
described by: 



r 



( 0.07 



g 
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)h 



<5.20) 



Equation <5.20) is based on empirical data. Kecman 
approximated the variation of the rolling radius along KA with 
a hyperbolic approximation: 



KA 




<5.21 ) 
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Where 1 is the distance -from K along KA . The linear 

k 

of the rolled length between K and A is given by: 


variation 


KA 

The energy along KA is: 


<5.22) 


2m z,KA 

W = — — 

KA 3r 

or 


<5.23) 


w=w +w "^fu2 2 2-1/2 

^7 V ^ ""la - 3 ^ ^A ’ 


<5.24) 


The relative rotation along yield lines PN, and 
whereas along KL, and MN the angle is given by: 


QM i s P , 


^A 

t, = arctan ( — ) 

^A 

Finally, the energy absorbed by these yield 

approximated by: 


lines is 


""e = ^ 2harctan(^)] 

' A 


<5.25) 


At this point, the total energy absorbed by the 
be expressed as: 


h i nge may 


8 

w(a) = Z w. (a) 
i=l ^ 


<5.26) 


This result is valid for a rotational angle less 
jamming angle. Kecman determines the moment-angle 


than the 
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relationship numerically -from the -following equation; 

M(a) = w(g + Aa) - VJ(g) (5.27) 

Ag 

Where Aa represents a small but -finite increment o-f the hinge 
angl e , a. 



This -formulation assumes that the material is ductile 
enough not to -fracture during hinge rotation or subsequent 
bending collapse. This is considered a sa-fe assumption -for 
common structural steels. The average plastic moment is given 
by the -following equation; 



m = 
P 




(5.28) 



The yield stress, and the maximum nominal -flow stress as 
-functions o-f temperature were determined in Chapter 5. 



5.1.5 The Final Form o-f the Hinge Moment 

The analysis in the preceding section is not valid -for 

the first phase of deformation. This phase will be 

approximated by assuming that for a moment greater than the 

fully plastic moment, M , is constant and equal to M . This 

o o 

would be a conservative assumption. For hinge angles greater 
than the jamming angle, Kecman offers the following formula; 

M(g) = M(gj) + 1.4 [M - M(gj)](g- a^) (5.29) 
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This formulation has been programmed in a M i cr osof t-BASI C 
program. The source code is given in Appendix B. Figure 5.5 
shows the results of the program. 



5^2— id^;k'duk.-S-^iici-LJl£CZbj.cJjLL-S-ld£.±hdd 

The primary difference between Kecman''s method and this 
method is that Hayduk and Wierzbicki considered not only 
i nex tens i onal deformation modes, but also extensional 
deformation patterns. The effect of these extensional modes 
will be fully developed in this section. 

Hayduk's and Ul i er zb i ck i •' s formulation results in a closed 
form analytical solution which predicts the mean crushing 
force. Reference 11 specifically addresses the problem of 
analyzing the crushing process of thin walled, plate formed, 
open structures with particular emphasis on L shapes and 
cruci forms. The I beam can be approximated by two T shapes or 
four L shapes. Lower and upper bound solutions for the mean 
crushing strength of cruciforms was obtained by considering 
various modes of deformation. According to Hayduk and 
Wierzbicki, extensional deformation accounts for at least one 
third of the dissipated energy. 
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The end products of this section will include the 
application of Hayduk"s and Wierzbicki's method to an I beam, 
and reviews of both references 11 and 12. 



5.2.1 Collapse Modes 

Hayduk and Wierzbicki considered out of plane 

deformations in which the deflections exceed 10-100 times the 

section thickness and are of the order of magnitude of the web 

height and the flange width. The process of collapse starts 

as the compression flange loses stability, either elastically 

or plastically at a critical load, P . Following the buckling 

c 

deformation, the flange can maintain added loads until a 

maximum load, P , is reached. Eventually, as the resistance 

max 

decreases, the structure collapses completely. The structures 
being examined here often regain strength after the initial 
loss of resistance. Subsequent load peaks reveal the 
development of additional folds. 
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5. 2. 1.1 Isometric Trans-f ormat i on Modes < I nex tens i onal 

De-format i ons) 

These modes are o-f the type used by Kecman to analyze the 
collapse mechanism. In this type o-f de-formation, the energy 
is concentrated over hinge lines, while the rest o-f the 
structure remains rigid. An assumption o-f i nex tens i onal 
de-formation is that material continuity is maintained. During 
the trans-f ormat i on o-f two sur-faces, the lengths o-f 
correspond! ng arc elements is preserved. Hayduk and 
Ulierzbicki show the progression o-f various de-formation paths. 
These are illustrated in Figure 5.6. 

In Figure 5.6 (a), and <b), the -final state is -formed 
i nex tens i bl y . This can only be done i -f the boundaries have 
su-fficient Freedom to deForm. The element in Figure 5.6 <c) 

is initially Folded along the stationary hinge line AB ; 
Finally the element is rotated extensibly, <reF. 11). 



5. 2. 1.2 Hinge Collapse Modes 

There are two types oF hinge collapse modes. The First 
is a propagating hinge line which occurs when the hinge line 
is Free to move with respect to material points, and the 
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material o-f one flange is effectively transmitted to another 
flange. Extensional deformation occurs in the vicinity of the 
hinge line. This mode is a char ac ter i st i c of thin walled 
structures consisting of two intersecting plates such as the I 
beam being considered in this study. 

The second hinge collapse mode is the formation of a 
stationary hinge line. This mode is char ac ter i zed by a hinge 
line fixed in the mater i al . Side f 1 anges- wh i ch are i n i t i al 1 y 
rectangular are transformed into trapazoidal elements. 
Therefore, extension occurs over an entire flange. This mode 
will not be considered in this analysis, <ref. 11). 



5. 2. 1.3 Assembled Collapse Modes 



Hayduk and Wierzbicki 


demonstrated the 


me thod 


Oi 


assembl 


i ng 


the folding patterns 


shown i n 


Figure 5 . 


6. The 


patterns 


are 


folded into more complicated 


symme trie, 


asymme tr i c 


and mi 


X 

a 


collapse modes which 


are shown 


in Figure 


5.7. 
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5.2.2 Anal ysi s 



For the purposes o-f this study, the analysis o-f a T 
section will be used; this analysis is considered to be 
equivalent to an I beam analysis. In re-ference 12, the short 
T cross section under consideration was subjected to 

compression with clamped boundary conditions. A model is 
presented in Figure 5.8. 

For the purpose o-f estimating the mean crushing -force, 

P , the T section is expanded into a plate and an angle as 
m 

shown in Figure 5.9. 

The plate will collapse through the bending zone limited 
by the two horizontal moving lines. The angle will follow the 
folding mechanisms described by Figure 5.10<b). The collapse 
modes in Figure 5.10<b) is consistent geome tr i cal 1 y and 
k i nemat i cal 1 y . The regions are described as follows: 

<I) Four plane trapazoidal elements as rigid bodies. 

<II) Two sections of cylindrical surfaces at which 
continuous bending takes place without any extension. 

(Ill) Two sections of conical surfaces in which material 
is bent and rebent again as the material moves from one 
flange to another. 
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Mccsl Of T 




Plate 




ArtQLE 



pOLPiHG, 

■piquRe £>. 10 



(Ik/') A section o-f toroidal sur-face which produces 



extension in a c i rcum-f ren t i al direction and continuously 
changes principal curvature in the other direction, (ref. 
1 1 ) 

The following assumptions will be used in this analysis; 

- Neglect the short side of the angle element. 

- The length of the two horizontal hinge lines is 
appr ox i mate 1 y equal to the length of the large side of 
the angle element. 

- Two incl ined hinge 1 ines are of a length equal to a 
half wave length, H. 

- The length of the toroidal surface in the vertical 
direction is neglected. 

The energy dissipated of an angle element is given by 
ummation of three terms, (ref 11): 

1. The energy of a continuous deformation field in the 
section of a toroidal shell; 

El =■■ 9.28M^Hr/t^ (5.30) 

2. The energy dissipated by the discontinuous velocity 
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field in the horizontal hinge line of length, 1; 



E 



2 



£ 

2 o 



<5.31 ) 



3. The work done by an inclined hinge line through the 
deformation process: 



E 



3 



2.22M H^/r 
o 



<5.32) 



where : 



r = The small radius of the toroidal surface. 

2 

M = (T’tt ] /4, the fully plastic yield moment per unit 
o o w 

width 

H = Mode half-wave length 



a — 0.7a 
o ou 

a = Ultimate strength 
ou 

t = Web thickness 



w 



For the models given here: 



2HP = E, + 4(E^) ^ + 4(E )^, + 2E 

m 1 2 web 2 flange 3 



<5.33) 



After substitution of the dissipated energies into <5.33): 



t. 2 



V^o = + d(:^) ] + 2,22^ 

w H „ 



<5.34) 
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Optimizing equation (5.34) by using the following partial 
derivatives: 



d m 
3r 



= 0 



<5. 35a) 



am 

3h 



= 0 



Equating (5.35a) and <5. 35b) yields: 



<5. 35b) 



r2.22 2, , ,^f 

r . I— (<= ^ d(— > 

4 • 64 V7 



TT^t 



2 



H - r W ^ ^ 1/3 

^i2.22) (4.64) ) 1 



w 



(5.36) 



(5.37) 



Ultimately, the three basic mechanisms of energy 
dissipation appear equally weighted in the final equation for 
mean crushing force: 



t 2 

P /M = 3 [ (4.64) (2,22)7r(l + k(-/-) ) (f-) ] 

mo t t 

w w 



(5.38) 



Equation (5.38) may be simplified as follows 



c m1/3 



P / M = 3[A,A^A, (— )] 
m o 1 2 3 t 

w 



(5.39) 



where : 



k = ^ 
c 



(5.40) 



= 4.64 



<5.41 ) 






^2 = ^[1 + k(— ) 



w 



A3 = 2.22 



(5.42) 



(5.43) 



Note also that equations (5.34), (5.36), and (5.37) may be 
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r 






« 



■ ^ 





simplified into the following form 



* ^37 

w 


<5.44) 


t(A A /A-^) (Ct 
Z J 1 w 


<5.45) 




<5.46) 



The application of these equations to the case of an I 
beam requires that the geometric parameters be modified. An I 
beam is two T sections, hence the dimensions are doubled; but, 
b = 2c, which yields the following relations. 

P /M = 3[A-A_A- (b/2t (5 47) 
1710 123 w 



k = 2a/b 



<5.48) 



It is now possible to apply this relation determined for 
the crushing case to the bending case. The geometry of this 
crushing mode is given in Figures 5.11 and 5.12. 



Recall the propagating hinge line mode as described in 
paragraph 5. 2. 1.2. This mode describes the global geometry of 
the crushing process and is determined purely by the width of 
the two side panels, d, which is prescribed, the wavelength, 
H, which is to be determined, and the crushing rotation angle, 
Q, which is a parameter of the process. The relation between 
the axial displacement, <^ , and the angle, is: 

6 = 2H(1 - cosB) <5.49) 
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Reference 12 determined an expression for the instantaneous 
crushing force, P, normalized with respect to P ; 



0. 9cosB 



T, j-n fl , ^ , U.ycosp 

m ~ 3 rrsinS tan6 (1 + sin^B) V 



Z-] 



<5.50) 



Figure 5.13 gives the geometry of the bending case. From 
Figure 5.13: 



* ^ 
tanj 



g/2 

b 



<5.51 ) 



Thus, it may be said that: 

P/P = f(6) = f(2btan^) 
m 2 



<5.52) 



For a section width, b: 



M = j b 



<5.53) 



M = ^ b 
m 2 



<5.54) 



Where M is the mean moment due to bending rotation, a, as 



m 



defined in Figure 5.14. After equating <5.53), and <5.54): 



: p 

m ni 



<5.55) 



Figure 5.15 gives P/P 



equations <5.49), and <5.50), 
equivalent description of 

<b/H) C tan<a/2) ] . 



as a function 


o+‘ 


S/2H based 


on 


< ref . 


12) . 


Th i s 


graph i s 


an 


M/M 


as 


a 


function 


o+ 



m 
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For this purposes of this study, an approximate relation 

•for M/M and a will be determined. Assume equation <5.56) is 
m 

a valid -function -for M/M : 

m 

M/M = d/(6/2H)" <5.56) 

m 

where ! 

D, n are constants 



Assume the new curve goes through points P , and P , in Figure 

1 2 



5.15. 



P, (6/2H, M/M ) = (0.30, 1.0) 
1 m 

P_(6/2H, M/M ) = (1.0, 0.24) 
^ m 

Solving for D , n ; 



<5.57) 

<5.58) 



D = 0.24 



<5. 5?) 



n = 1 . 2 

Assuming n = 1.0, the following equation results: 



<5.60) 



F i gur e 
exact 
<5.61 ) 



" 6/2H 



5.16 shows how the approximate 
solution. Substituting equation 
, the following results: 



<5.61 ) 

solution compares to the 
<5.51) into equation 
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.ZC> -IS to 

3 , 

COHPAPJSOH OP ApPR.OXlMATe' 
Ap^ EcACT SiCuuTlO'K^ 



-pltSiuRe S', \ 6 



<5.62> 



M/iM = 



0.24 



m (b/H) tan (a/2) 
Considering equation (5.54); 



P b 



2 (b/H) tan (a/2) 



-] 



<5.63) 



Now substitute equation (5.47) into (5.63) 



M - l3HJA^A2A3(b/2t„)//^nb/21[2j|i^l 



(5.64) 



But M is the fully plastic yield moment for the flange 
o 

2 

a t 

M = — 

o 4 



(5.65) 



where t = t = t . 

f w 



The desired final result of this analysis is a relation 

between M and M . Rewriting equation (5.65): 

o 

M * = G / zdA 
o o A 

Recall equations (2.2) and (2.4): 

= -20. 5T f 40,000.0 
2 

zdA = abt + ^ 

Finally, substitute equation (2.4), (5.66), (5.65), and 

(5.46) into (5.64); 



(5.66) 

( 2 . 2 ) 

(2.4) 
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M/M ' = 
o 



3a t 
o 



4a^(abt + b t/4) 



0.24[(A2/A A )b^2t]^/^ 
'^^^2^ ^ f btan(a/2) ^ 



The simpli-fied -final result of this relation is 



, ^ 0.36TTt(b -I- 2a) (5.68) 

tan(a/2) (b(b 4a) ] 

Figure 5.17 shows graphically the resulting momen t-angl e 
rotation relationship as determined by Hayduk and Wiersbicki. 
Appendix C contains the source code for a Microsoft BASIC 
program which determines the moments resulting from equation 
<5.72> . 



2^3— £omp^J.£.aB— £ja^it2±^_jo£— JLb£— Jwo— tlelbcuds 

The primary difference between the two methods of 
analyses is that Kecman considered only i nex tens i ona 1 
deformation modes, whereas Hayduk and Wierzbicki considered 
extensional deformations as well as i nex tens i onal 
deformation. Since the second method considered more energy 
dissipation than the first, the allowable moments due to 
deformation are greater. However, for angles greater than one 
tenth of the jamming angle, Kecman ''s method yields greater 
moments. The errors induced by his numerical integration 
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I 



1 



technique are compounded as the angles increase. Figure 5.18 
shows graphically how the results of each analysis differs. 
The dimensions of the section analyzed are given below: 

a = 12 inches 

b = 12 inches 

t = t = .6 inches 
f w 

The maximum moment given in Figure 5.18 is the fully plastic 

yield moment of the section, hereafter called M , which is 

o 

only a function of geometry and temperature of exposure of the 

beam. Hence, M is not a function of the analysis method. 

o 

Clearly, as the angle of rotation increases the moment 
decreases; and as the temperature increases, the moment 
decreases. Appendices B and C also contain results which show 
the effects of changing the geometry of the section. For both 
analyses, the moments are decreased as the section thickness 
decreases; the web height decreases; or the flange width 
decreases . 

These methods of analyses are considered accurate 
predictions of the failure mechanisms of a beam exposed to 
fire. Recommendations for using the information provided by 
these analyses will be given in the following chapter. 
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Chapter 6 



CnNCl.llSIONS AND RECOMMEtiDftllDt^ 



A summary o-f each chapter is presented here with 
concluding remarks. Finally, recommendations are of-fered -for 
using this analysis -for -future studies. 



i£_JL— Summajix 

As stated in the introduction, the ultimate objectiue o-f 
this study was to show how much the strength o-f a shipboard 
structural member is diminished when exposed to fire. To this 
end, Chapters 2, 3, and 4 provided input into the analyses of 
the failure mechanisms studied in Chapter 5. The results of 
the preceding chapters will be summarized here. 

Chapter 2 developed the model of the structure. Simple 
observations revealed that the ship is a redundant structure 
composed of beams, plates, and columns. The motivation 
towards choosing an I beam for the analysis was developed. 
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The fTiciterial selected -for the model was ASTM-A36 steel as its 



physical and mechanical properties are similar to the 
properties o-f ABS Grade A steel, a common shipbuilding steel. 
The equations -for the temperature dependent mechanical 
properties were the output o-f this chapter to be used in the 
analysis. 

Chapter 3 provided background -for char ac ter i z i ng 
shipboard -fires. Three parameters which are assumed to be 
accurate predictions o-f the contributing -factors of a fully 
developed fire were defined. These factors are; the duration 
of the fully developed firej the overall penetration flux; and 
the average fire temperature. Since fire is indeed a complex 
process with several unquan t i f i abl e variables, a time 
temperature curve was used as a definition of the fire 
process. For the purposes of this study, a step function to 
o 

1700 F was used as the time temperature relation. 

Chapter 4 was the heat transfer analysis which determined 
the relationship between the elevated temperature exposure and 
the distance along the length of the beam. The results of 
this analysis led to the conclusion that a beam exposed to 
fire will experience the same failure mechanisms as a beam 
under a point load. 

Finally, the primary thrust of this study was to analyze 
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the -failure mechanisms due to fire exposure. Two methods were 
fully developed in the preceding chapter and the results were 
summarized graphically in Figure 5.18. 



Eacommanda±xons 

The ultimate objective of any structural design process 
is to develop a structure which is strong enough to withstand 
all expected service loads and a high percentage of unexpected 
loads. The information provided by the analysis in Chapter 5 
is considered useful input for the design of shipboard 
structures. It is possible through the use of Figure 5.18 to 
determine how the stiffness of a structural member is 
decreased when exposed to elevated temperatures . The 
diminished stiffnesses may be used as input to a finite 
element program which simulates a ship structure. The effect 
of a shipboard fire upon the entire structure may then be 
de term i ned . 

The examination of the effect that one weakened member 
has on the other members of a redundant structure is 
considered fertile ground for future studies. Reference 16 
may provide insight into the behavior of this phenomenon. 
Zubelewicz and Mroz numerically simulated rock burst processes 
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treated as problems of dynamic instability. According to 
reference 16 , local plastic failure occurs at one location 
(such as in a member of a Shipp's structure) and the elastic 
energy flows to the hinge and is the driving force into the 
damage. Energy then flows into the neighboring hinges, which 
is to say that one unstable member may produce a chain 
reaction and other members may subsequently fail. 

As a final word, it is highly recommended to utilize the 
information provided in this study as an aid for examining the 
effect of fires on the ship^s structure. Ultimately a ship 
may be designed to withstand shipboard fires. 
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Appendix A 



tl£AI_I£ANS£££_ANAL:i:SiS_£EDGSAtl 



10 R£M PROGRAM TO CALCULATE TEMPERATURE GRADIENT IN A BEAM 
20 DIM TEMP<20) ,Y<20) ,2<20) ,TOP<20) 

25 INPUT "ENTER THE BEAM LENGTH IN FEET" ;XL 
30 INPUT "IS BEAM SOLID <S) OR I <I)";N$ 

40 IF N«="I" THEN GOTO 90 
45 LPRINT " SOLID BEAM" 

50 INPUT "ENTER WIDTH OF BEAM" ;W 
60 INPUT "ENTER DEPTH OF BEAM";D 
65 LPRINT " WIDTH, DEPTH=" ; W;D 
70 B=SGR< < 12/5.4)#< < l/D)+< 1/W) ) ) 

80 GOTO 155 
84 LPRINT " I BEAM" 

90 INPUT "ENTER WIDTH OF FLANGE" ;W 
100 INPUT "ENTER DEPTH OF WEB";D 

105 LPRINT " WIDTH OF FLANGE, DEPTH OF WEB=" ; W;D 
110 INPUT "ENTER THICKNESS OF FLANGE" ;TI 
120 INPUT "ENTER THICKNESS OF WEB" ;T2 

125 LPRINT " THICKNESS OF FLANGE, THICKNESS OF WEB=" ;T1 ;T2 
130 P=2*<D+2*W+2*T1-T2) 

140 S=2*<T1 *W)+<D*T2) 

150 B=SQR< < 12/5.4)*(P/S) ) 

155 LPRINT " LENGTH OF BEAM=";XL;" FEET" 

156 LPRINT " B=";B : LPRINT " " 
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160 BL=B*XL 



170 BOT=EXP<BL)-EXP(-BL) 

ISO X=XL/20 

IPO LPRINT " PT NO X T-Ta/Tb-Ta" 
200 FOR KX=1 TO 20 
210 Y<KX)=X»KX 
220 2<KX)=B*<XL-Y<KX)) 

230 TOP(KX)=EXP<2<KX) )-EXP<-2<KX)) 
240 TEMP<KX)=TOP<KX)/BOT 
250 LPRINT " " ;KX,Y<KX) ,TEMP<KX) 
260 NEXT KX 
270 END 
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Appendix B 



K£CtlAfcJ^S_ANAJ_:»rSlS_Pfi£lE£Ald 



10 DIM W<20) 

20 INPUT "TEMPERATURE OF EXPOSURE, DEGREES F" ;T 
30 INPUT "FLANGE WIDTH, INCHES" , -A 

40 INPUT “WEB HEIGHT, INCHES" ;B 

41 INPUT "THICKNESS, INCHES" ;TH 

42 LPRINT " " 

50 LPRINT " " 

51 LPRINT " TEMPERATURE=" ;T 

52 LPRINT " FLANGE WIDTH=";A 

53 LPRINT " WEB HEIGHT=";B 

54 LPRINT " THICKNESS=" jTH 
60 S0=-20 .5*T+40000 .0 

65 IF T>=1?51 THEN S0=4 . 5 

67 M0= S0*B*TH#<A+B/4) 

68 LPRINT " FULLY PLASTIC MOMENT, INCH-POUNDS" ;M0 
70 J=0 

90 U=3. 14159265 
100 IF A <= B THEN H=A/2 
110 I F A > B THEN H=B/2 
120 Q1=(H-.5*TH)/B 

130 D=2*<ATN(Q1/SQR(-Q1*QH-1))) 

131 LPRINT “ JAMMING ANGLE, RADIANS" ;D 

132 LPRINT " ANGLE; MOMENT/MO; MOMENT" 
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140 X=SU*<TH''2)/4 

150 FOR K=D/10 TO 11*0/10 STEP D/10 

151 PI = K/2-K''3/48+K''5/3840 

152 P2 = l-K”' 2/8+ K'' 4/384 

153 P3 = B*P1 

154 P4 = (2*H - P3) 

155 P5 = B*P2 

156 P6 = H*P1 

157 P7 = P1/P2 

165 YB = P5-<SQR<P3*P4) ) 

170 R=( .0 7-h0''70)*H 

180 2A=<P3*P1 )-P6+(SQR<P3*P4) >*P2 

1?0 Q1=1-<B/H)*P1 

200 Q2=ATN<Q1/SQR<-Q1*Q1 + 1 ) ) 

210 Wl=2*A*<U/2-K/2-Q2) 

220 W2=A*<U-2*Q2) 

230 W3=( P3*P1 -P6+ < SQR( P3*P4) ) *P2) *( U-2*Q2) 
240 W4=H*U 

250 YC=<H*P7+P5-<SQR<P3*P4) ) )/<l+P7''2) 

260 XC=YC*P7 

270 Q1=2A/<SQR< < ( H-XC) ''2) + < ( YC-YB) •'2) ) ) 

230 Q3=ATN(Q1) 

2?0 W5=2*B*Q3 

300 W7=4/3*<2A/R)*<SQR< ( H''2) + < YB''2) + ( 2A''2) ) ) 
310 W6=2*(H/R)*2A 
320 Q1=ZA/YB 
330 Q3=ATN<Q1) 
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} 



II 



340 W8=2*<A*<K/2)+2*H*Q3) 



350 J=J+1 

360 UK J)=U)1 +UJ2+W3+UI4+UI5+UI6+W7+UI8 

370 NEXT K 

380 FOR IX=1 TO J-1 

390 Y=<TH/<B*<4*A+B) ) )*<W< I'/.+ l )-W< IX) )/(D/10) 

391 M1=Y*M0 

400 LPRINT ■' " ; IX*D/10 ;Y;M1 

401 NEXT IX 

402 LPRINT ” '• 

403 LPRINT '■ " 

404 LPRINT ■■ " 

410 LPRINT '• “ 

420 END 



103 



Appendix C 



iJAYDJJK_ALJD_yi£RZBiCKJ.^£_AiiJALYSlS_£E£IE£Ad 



10 INPUT "TEMPERATURE OF EXPOSURE, DEGREES F" ;T 
30 INPUT "FLANGE WIDTH, INCHES" ;A 
40 INPUT "WEB HEIGHT, INCHES" ;B 
50 INPUT "THICKNESS, INCHES" ;TH 

61 LPRINT " " 

62 LPRINT " " 

70 LPRINT " TEMPERATURE " ;T 
SO LPRINT " FLANGE WIDTH" ;A 
90 LPRINT " WEB HEIGHT" ;B 
100 LPRINT " THICKNESS" ;TH 
120 Tl=<A*B#TH+B*B*TH/4) 

130 T2=-20 .5*T+40000 

135 IF T>=1951 THEN T2=4.5 

137 M0=T2*T1 

138 LPRINT " FULLY PLASTIC MOMENT, INCH-POUNDS" ;M0 
140 T3=0 .36*3. 141 593*TH*<B+2*A) 

150 T4=1/<B*(B+4*A) ) 

170 LPRINT" ANGLE; MOMENT/MO; MOMENT" 

171 IF A<=B THEN J=A/2 

172 IF A>B THEN J=B/2 

175 Ql=< J-.5*TH)/B 

176 D=2*<ATN<Q1/SQR<-Q1*Q1+1 ) ) ) 

180 FOR K=D/10 TO D STEP D/10 
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181 T5=TAN<K/2) 



190 M=T3*T4/T5 
195 M1=M*M0 

200 LPRINT " '';K 

201 NEXT K 

202 LPRINT " ■■ 

203 LPRINT " " 

204 LPRINT " " 

210 LPRINT " " 

220 END 
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